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Abstract 

We consider the fc-folded sine-Gordon model, obtained from the usual version by 
identifying the scalar field after k periods of the cosine potential. We examine (1) 
the ground state energy split, (2) the lowest lying multi-particle state spectrum and 
(3) vacuum expectation values of local fields in finite spatial volume, combining the 
Truncated Conformal Space Approach, the method of the Destri-de Vega nonlinear 
integral equation (NLIE) and semiclassical instanton calculations. We show that 
the predictions of all these different methods are consistent with each other and in 
particular provide further support for the NLIE method in the presence of a twist 
parameter. It turns out that the model provides an optimal laboratory for examining 
instanton contributions beyond the dilute instanton gas approximation. We also 
provide evidence for the exact formula for the vacuum expectation values conjectured 
by Lukyanov and Zamolodchikov. 
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1 Introduction 

Completely integrable 2 dimensional quantum field theories are important sources of non- 
perturbative information but are rather rare. Thus there is a considerable interest in 
modifications or deformations that generate new integrable models from an old one. The 
extensively studied examples of such modifications include the integrable deformations of 
CFTs III and the class of boundary integrable theories where the new boundary condi- 
tions preserve integrability. 

In this paper we consider a slightly different integrable modification of sine-Gordon 
(SG) model, which was first mentioned in Q and discussed in some detail in 0. In this 
model the period of the sine-Gordon field (f) is k times the period of the potential for some 
generic k e NQ. Thus the new model — which we denote by SG(/3, k) — contains the folding 
number as a parameter in addition to the usual coupling constant p. All properties of 
the ordinary SG model (SG(/3, 1)) that rely only on the local properties of the field 0, like 
classical integrability, or the existence of conserved higher spin quantities, will also hold 
for SG(/5, k). Thus we have a rather clear intuitive picture of the various (infinite volume) 
excitations of SG(/9, fc), at least for < Svr. This picture predicts the exact S'-matrices 
of the scattering among the particles corresponding to these excitations in terms of the 
well known S'-matrix of SG. Nevertheless the spectrum of SG(/5, k) is not identical to that 
of the ordinary SG, refiecting the importance of the boundary conditions imposed. In 
particular the quantum theory SG(/5, k) has a fc-fold degenerate vacuum corresponding 
to the 'unidentified' minima of the potential. As a consequence SG{P,k) contains kinks, 
i.e. particles that interpolate between different vacua and have nontrivial restrictions on 
their multi-particle Hilbert space. Since these restrictions are k dependent they do give 
rise to differences between the theories with different k, that can manifest themselves in 
their finite volume spectra. 

In this paper we investigate SG(/3, k) in finite volume. We study three sets of problems 
in some detail, namely the split in the vacuum energy levels, the spectrum of the low lying 
multi-particle states, and finally the vacuum expectation values of exponential fields. In all 
three cases we compare the theoretical predictions with numerical data obtained by using 
the Truncated Conformal Space Approach (TCSA) 0. 

We derive the split of the vacuum energy levels by two methods, on the one hand we 
obtain it from the nonlinear integral equation (NLIE) appropriately generalized 
to describe SG{P,k), while on the other we perform an instanton calculation. These two 
methods give identical results for the leading part of the split, and the TCSA data show 
an excellent agreement with this prediction. This agreement confirms the correctness of 
the generalized NLIE. We also show that the TCSA data make it possible to extract and 
analyze the nonleading part of the split. 

We determine the volume dependence of the energy levels of the low lying multi-particle 
states by using the formalism of [||. This method relies heavily on the conjectured S- 

^This procedure can be carried out in any scalar field theory — not necessarily integrable — in which the 
scalar potential is a periodic function. 
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matrices describing the mutual scatterings in the multi-particle states. For large volumes 
this simpler and approximate method can be shown to give identical results to the exact 
NLIE, and we chose it because it has a clear interpretation in terms of the spectrum of the 
model. The — k dependent — degeneracies and the volume dependence of the multi-particle 
energy levels obtained this way agree very well with the TCSA data, and this agreement 
confirms that the conjectured 5'-matrices are indeed the correct ones. 

The study of the vacuum expectation values of exponential fields is motivated by the 
fact that in an explicit formula was proposed for this quantity at least in the ordinary 
sine-Gordon theory in infinite volume. We translate this explicit expression into the fi- 
nite volume SG(/3, fc), and compare it to the expectation values measured in the ground 
states found by the TCSA. The TCSA data are good enough to distinguish between the 
semiclassical and exact expressions given in 0, favoring the latter. The agreement we find 
proves two things. On the one hand it confirms the expression given in on the other 
it shows that the expectation, that everything, which in sine-Gordon theory follows only 
from the local properties of the scalar field and the Lagrangian, remains true also in the 
/c-folded model, is indeed correct. 

The paper is organized as follows: in Section 2 we describe the Lagrangian, the sym- 
metries and the local operators of SG(/3, /c), and recall the basic properties of TCSA as 
applied to this model. Section 3 contains the generalization of the NLIE. We investigate 
the vacuum structure of SG(/3, k) in finite volume in Section 4. Section 5 is devoted to the 
study of the multi-particle energy levels, and we analyze the vacuum expectation values of 
exponential fields in Section 6. We make our conclusions in Section 7. The paper is closed 
by an appendix where we describe in detail the computation of a determinant needed to 
complete the instanton calculation of the split between the vacuum energy levels. 

2 The /c-folded sine-Gordon model 

2.1 SG{(3,k): Lagrangian and symmetries 

The action of sine-Gordon theory in a finite spatial volume L is 

For later convenience, we also define a new parameter p with 



P 



To define the fc-folded theory SG(/3, k) we take the sine-Gordon field ip as an angular 
variable with the period 

^^^ + —k. (2.2) 
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This implies the following quasi-periodic boundary condition for the field: 

ip{x + L,t) = ip{x,t) + —km , m & Ia . (2.3) 

P 

The classical ground states are easily obtained: 

'^n = ^n, n = 0,...,k-l (2.4) 



which shows that the condition ( |2.2| ) corresponds to identifying the minima of the cosine 
potential with a period k (see Fig. ||). All classical solutions of the ordinary SG model 
are also solutions of SG(/3, /c), as the equations of motion are identical. However the 
static soliton solution, which, in the pure SG framework, on account of the identification 
(f = (f + '2,71/(3, interpolates between the same minimum, is now connecting different, 
neighbouring minima; i.e. the SG soliton becomes a kink in SG(/?, k). 

In the infinite volume (L = oo) quantum theory these ipn correspond to the vacuum 
states \n) which have the property 

{n\Lp{x,t)\n) = '^n . (2.5) 
P 

These states are all degenerate in the classical theory and also at quantum level when 
L = oo; however, tunnelling lifts the degeneracy in finite volume L < oo. 

Let us now examine the relevant symmetries of the action. One can define a Z^. group 
action generated by a unitary operator T in the following way: 

Tif{x,t)T-' = ipix,t)~j. (2.6) 

Then we have 

T\n) = \n + l mod k) (2.7) 

and the Hamiltonian 

H = j'^'^dx Qvr^ + 1 {d^^f + ^ (1 - cos/3^)^ , 7r(x,t) = d^vM (2.8) 
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commutes with T. The eigenvectors of T are the "Bloch waves": 

1^-) = 7r E 1^) ' = ' (2-9) 

^ m=0 

T|^9„) = e-^''"|^?„) , (2.10) 

which are eigenstates of the Hamiltonian as well. In fact, these states can be continued to 
finite volume and are eigenstates of if when L < oo, while there are no natural counterparts 
of the states \n): for finite L we define them using the inverse formula 



|m) = ^ 5^ e— "'^ . 



fe-i 

One can also introduce a Z2 transformation S which is defined as 

Sifix, t)S-^ = -(^(x, t), = S-^ = S. (2.11) 

It acts on the ground states as 

S \n) = \k-n) = \-n) , S = = \^k-i) , (2.12) 

and commutes with the Hamiltonian H. The S and T transformations together generate 
the discrete group D^.. 



2.2 The spectrum of local operators 

As a consequence of (|2.2| , |273|) , the exponential fields^ 

= exp (iPj^) (2.13) 
are well-defined local operators provided m G Z. One can easily compute 

TVmT~' = e-'^^Kn , SVmS-' = V.m . (2.14) 

In the short distance limit, the behaviour of the correlation functions is described by a 
c = 1 compactified free boson with the Lagrangian density 

C=^d^xd''x- (2.15) 

OTT 

^In writing Vm we assume that the exponential fields are normal ordered and are normaHzed such that 
the short distance limit of their two point function is 

(V„i{x)V^,niy)) \x - y\ k-M^ for \x - y\ ^ . 
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In order to determine the complete spectrum of local operators, let us consider this UV 
limiting theory. The field x lives on a circle with compactification radius r 

X ~ X + 27™r , n G Z 



where r can be determined from (|2.2|) . Taking into account that the normalizations of 
and X differ by a factor of v^ivr we get 



47r , , , 

— k. (2.16) 

The above theory has a U{l)j^ x f/(l)^ Kac-Moody symmetry. The primary fields under 
this algebra are vertex operators V(n,w) with left/right conformal weights given by 

_i_ 1 /n , wr\" 



Ar„,„,= 2(-±^j . (2.17) 

n is the field momentum quantum number, while w is the so-called winding number. The 
local operators Vm correspond to V(m,o)- The sine-Gordon potential can be identified as 

:cos/5(^:= ^(V(fc,o) + V^-kfi)) ■ (2.18) 

There are only two possible maximal local operator algebras in a c = 1 free boson theory 
with compactification radius r f^: 

A = {V(„,^) : w e Z, n e Z} , 

Af = {V(n,w) ■■ w eZ, neZ + w/2} . (2.19) 

The first one corresponds to a bosonic model, while for the second one the operators 
corresponding to w odd are fermions. This gives the complete list of possible local theories 
in the "sine-Gordon class" where standard sine-Gordon/massive Thirring corresponds to 
choosing k = 1 and the algebra Ab/Af , respectively. 

2.3 Truncated Conformal Space for SG(/?, k) 

In order to support our theoretical considerations we shall use the well-known Truncated 
Conformal Space Approach (TCSA) to obtain numerical data for the model SG(/3, fc). 
This method was developed by Yurov and Zamolodchikov ||5| for perturbations of Virasoro 
minimal models and more recently extended to perturbations of c = 1 theories in ||10| 
where we refer the interested reader for more detail. 

Here we limit ourselves to recalling the basic facts. We can represent the Hamiltonian 
as an infinite hermitian matrix on the space of states of the c = 1 conformal field 
theory built on the algebra Ab ( |2.19|) P|: 



n= (2.20) 



^From now on we will restrict ourselves to the model based on the algebra At, as its fermionic coun- 
terpart is very similar. 



2 THE K-FOLDED SINE-GORDON MODEL 



7 



where 

^(n,m) = span {a_fci • • • \V(n,m)) , ki,...,ki E Z+} (2.21) 

is the Fock space built over the primary state |V(„,m)) with the negative frequency modes 
of the conformal free boson field x (|2.15| ). The Hamiltonian takes the form 




H = ^ \ Lo + Lo - ^ld + X^^B \ , (2.22) 

where Lq and Lq are diagonal matrices with their diagonal elements being the left and 
right conformal weights, Id is the identity matrix, 

h = ^ = — ^ 2.23 

An p+1 ^ ' 

is the scaling dimension of the perturbing potential and the matrix elements of B between 
two states |$) and are 

= ^ ($1 V(fc,o)(l, 1) + V(_fc,o)(l, 1) 1^) . (2.24) 



The parameter A is connected to Hq in eqn. (|2.1|) via a known relation ||ri|. The matrix 
elements of B can be calculated in closed form. We would like to call attention to the 
fact that the matrix elements of all vertex operators V(n,m) in the basis (|]2l]) are real 



numbers. As a result the Hamiltonian is a real symmetric matrix which will be important 
in establishing (|6^) . 

We choose our units in terms of the soliton mass M which is related to the coupling 



constant A by the mass gap formula obtained from TEA in |Tl|: 

A = K{h)M'^-^, (2.25) 

where 

(Note that we use the same massgap relation in the fc-folded model as in ordinary sine- 
Gordon, due to our previous argument about the relation between local properties in 
the two models). In what follows we normalize the energy scale by taking M = 1 and 
denote the dimensionless volume ML by /. For numerical computations, we shall use the 
dimensionless Hamiltonian 
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We diagonalize the matrix in a truncated Hilbert space defined as 

HTCsis, W, Ecut) = { 1^) : (^0 -Lo)\^) = s\^),Qk\'^)=w\^), HcfT |^) < ^cut 1^)} , 

(2.28) 

where, besides imposing an upper bound on conformal energy, we restricted the Hilbert 
space to a given value of the conformal spin Lq — Lq and of the topological charge Qk (the 
winding number of the free boson x) since these commute with the Hamiltonian H. As 
in the repulsive regime (p > 1) the TCSA for the SG theory is plagued by UV problems 
10|| , in this paper we restrict all of our numerical studies to the attractive regime (p < 1). 



For the purposes of producing the numerical data used later in the paper we typically used 
values of ii^cut that correspond to roughly 6000-12000 states. 



3 SG(/?, k) in the NLIE framework 

In a finite spatial volume the spectrum (the ground and the excited states) of sine-Gordon 
theory ( |2.1j) is described by a nonlinear integral equation (NLIE) P, 0. We shall start 
with a more general equation than in the normal sine-Gordon situation by introducing a 
twist angle ^ a \a Zamolodchikov [0], originally motivated by considerations related to 
polymers. Later it appeared in the description of the finite volume spectrum of Virasoro 
minimal models perturbed by $(1,3) ||13|- It corresponds to switching on a chemical potential 




coupled to the topological charge ||14|. The full twisted NLIE reads 

/oo 
dx G{X~x- ir]) log (1 + (_i)^e^^(^+*'^)) 
-00 

/oo 
dx GiX-x + 17]) log (1 + (_i)5e-'^("-*'')) (3,1) 
■00 

where M is the soliton mass, L is the volume, r/ is a suitably chosen real shift, 6 is either 
or 1, and the kernel G reads 

sinh ^(P-^)^ 

dfce^^^ ^ r. (3.2) 

2 sinh ^ cosh f ^ ' 

The twist angle d can be restricted to lie in the range — vr < t9 < vr without loss of generality: 
this choice simplifies the description of the results of the UV calculations. 

The term g{\\\j) is the so-called source term, composed of contributions from the holes, 
special objects (roots/holes) and complex roots. The form of g is specific to the state in 
the spectrum one wants to describe; for states with no particles, = (at least for L large 
enough, where the so-called special objects do not appear). We denote their positions by 
the general symbol {Aj} = {/i^ , yk , Ck , Wk} (h stands for holes, y for special objects and 
c (w) for close (wide) complex roots). The source term takes the general form 

Nh Ns Mc Mw 

9{M^j) = X] X(A - /ifc) - 2 ^ X(A - 2/fc) - X] ~ ^*:) ~ 5Z ~ ^'^^ii ' 

k=l k=l k=l k=l 
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where 

-A 



X(A) = 27r / dxG{x) (3.3) 
Jo 

and the second determination for any function /(A) is defined by 



/l^Jii I /(A) -/(A-msign(53mA)p) , p<l, ^ ' 

whenever |53mA| > min(7r,7rp). 

The source positions A^ are determined from the Bethe quantization conditions 

Z{\,)=2nl,, I^eZ + ^, (3.5) 

where Ij are the Bethe quantum numbers (for wide roots one must use the second deter- 
mination of Z, defined as in ( p.4| )). 

Given a solution for Z, the energy and momentum of the state can be computed from 
the formulae 

Nh Ns 

E = Eh^i], + M ^ cosh hj - 2M cosh yj 

Mc Mw 

—M coshcj — M ^^(coshwj)n 



S iZ(x+irj) 



/oo dx 
— 23m rsinh(x + irj) log(l + {-ife 
-oo 2vr 



Nh Ns 

M sinh hj - 2M ^ sinh yj 

j=i j=i 

Mc Mw 

—M sinh Cj — M ^^(sinh Wj)ii 

/°° rjr 
—253m rcosh(x + 17]) log(l + 
-oo 2vr ^ 



^buik = - Jm^L tan ^. (3.6) 



where 



For more details on the NLIE for excited states we refer to the literature [10|, |13|, |15|, |16| . 

A detailed calculation of the UV limit of the twisted NLIE ( |3.1|) was performed in 
T3|| . For our case, remembering that the identification of the perturbing potential (|2.18|) is 
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different from that of Feverati et al., we have to trace the appearance of the sine-Gordon 
twist parameter k in their formulas. Provided one chooses 



6 = Nh mod 2 



(3.7) 



the ultraviolet limit is described by states in the conformal family of vertex operators of 
the form 



V, 



{n,m) 



nekZ+ — 
2tc 



m=-{NH- 2Ns - Mc - 2Mwe{p - 1)) , 



where 9{x) denotes the Heaviside step function. 

Therefore to describe the UV spectrum Ab (|2.19|) of the fc-folded sine-Gordon theory 
one must choose the following values for the twist parameter: 



27Tn 






'k' 




[-H 










, . . . , 


2 



(3.8) 



which is mod 27i equivalent to the set of values in (|2.9|) . 
We also have to restrict 



Q = Nh- 2Ns -Mc- 2MwO{p - 1) e kZ . 



(3.9) 



The quantity Q counts the topological charge in the units of the usual sine-Gordon theory 
{SG{P,k = 1)). The formula (|3.9|) expresses the fact that all states of SG{P,k) which 
satisfy the periodic boundary conditions must be compatible with (|2?3|) as a consequence. 



4 Vacuum structure and instantons in finite volume 

In this section we investigate the vacuum structure of the fc-folded sine-Gordon model 
and show that it provides a laboratory to further test the NLIE as well as for analyzing 
the higher order corrections to the dilute instanton gas approximation (DIG). We test the 
NLIE by comparing its predictions both to the results of the instanton calculus and to the 
TCSA data. The second possibility is due to the fact that terms beyond the DIG give 
contributions to certain quantities, where they are not suppressed by lower order terms. 
First we give a group theoretical, thus qualitative description of the vacuum spectrum, then 
we analyse it quantitatively using the NLIE and finally we make the instanton calculation. 

4.1 Symmetry considerations 

In analyzing the vacuum structure we start with the L — > oo limit. In this case Derrick's 
theorem forbids the existence of instantons. We have k different degenerate minima with 
the corresponding states given by |m), m = 0, . . . , /c — 1. The symmetry transformations 
act on |m) by eqns. (|2.71) , ( p.l2| ) and they commute with the diagonal Hamiltonian. 
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Now we consider the theory for finite L. For finite L instantons exist and they lift 
the degeneracy of the ground states. As a consequence we have a unique vacuum state 
invariant under D^. and the Hamiltonian is no longer diagonal in the basis spanned by |m). 
Nevertheless from its symmetry properties we can determine its form. Since it commutes 
with the symmetry transformations it must belong to the center of the group algebra of 
Dfc. The generators of the center can be obtained by summing up all elements of a given 
conjugacy class with the same weights. One can show by direct computation that in the 
particular representation considered the elements of the center which correspond to the 
S transformation can be expanded in terms of those generators which correspond purely 
to the T transformation. Consequently, the most general element of the center has the 
following form: 

H = Eo{L)ld + Ei{L){T + T-') + ... + Ei{L){T + T~') + . . . . 

We can diagonalize this matrix by diagonalizing its commutant, i.e. the representation 
matrices of T itself. The eigenvectors are the states \'dm) introduced in eqn. and the 
corresponding eigenvalues of the Hamiltonian are: 

= Eo{L) + Ej{L) cos ( —jm j , m = 0, . . . , A; - 1 , 2Ej = Ej . (4.1) 
i=i ^ ^ 

Thus for finite L, instead of the fc-fold degenerate ground states we have the following 
vacuum structure: there is a nondegenerate lowest eigenvalue, Hq, while the rest of the 
eigenvalues come in pairs, Hr = Hk_r for 1 < r < [A;/2], at least for odd k. For k even 
Hkj2 is non-degenerate as well. Clearly from the knowledge of the energy levels Hm we can 
recover all the coefficients Ej{L). 



4.2 Leading finite size corrections to the vacuum energy 

The NLIE for the twisted vacuum reads 

/oo 
dx G{\-x- if]) log (1 + e^^(^+^'^)) 
-oo 

/oo 
dxG(A-x + ir/)log(l + e-^^(^~*'')) , (4.2) 
-oo 

Here we relax the condition ( p.8| ) on in order to keep the calculation more general. As 
7? + 27r is a symmetry of the NLIE ([1.2|) we can restrict the value of the twist angle {} 
to — vr < {} < TT. The value = corresponds to the untwisted sector. The vacuum energy 
can be obtained using 

/°° dr 
— sinh(x + zr/)log(l + e*^(^+^'')) , (4.3) 
■oo 27r 
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where the counting function Z{x) is a solution of the vacuum NLIE ([1.2| ) and we omitted 
the universal bulk energy term ( p.6|) . For large values of the volume L, the solution is 



Z{x) ^^ + 1 smh{x) , l = ML . (4.4) 

The value of rj must lie inside the analyticity strip for the kernel G, i.e. \ri\ < min(7r, 7rp) 
("first determination"). However, there is no singularity whatsoever in the integrand of 
E{L) when using the above approximation for Z and the contour can be shifted to rj = m (2 
for any value of p (in fact, the result is independent of p itself). One obtains 

— cosh(x) [log (1 + 6^''-'^°^^") + log (1 + e-^^^-'^^-hx^j _ ^4 5^ 

Expanding the term log(l + . . .) in Taylor series, a short computation gives 

2A/f °° f— 1')"^"'^ 

E(^,L) = KUnl) cos n^, (4.6) 

71 ^-^ n 

n=l 

where Kx{z) is a modified Bessel function of the second kind: 

/"°° da; 
= \ — cosh(a;)e 



- z cosh X 



We remark that this result is exact at the free fermion point p = 1 (where G = due to 
( ^^ ) and it can be shown to resum into 



a / /2 \ 

-5^K/((2m-l)7r + ^)2 + /2-(2m-l)vr-— — + 

TT ^-^ \ l\lm — 1 TT / 

m=l ^ \ J / 

\ t (v/((2™ - 1). - + P - (2ra - 1). - ^p^) } (4.7) 



with C = 0.57721566 . . . being the Euler-Mascheroni constant. This is exactly the result 
for a free Dirac fermion with twisted boundary conditions in finite volume L. 

4.3 Instantons in finite volume 

With the exception of the leading n = 1 term in ([4.6|) , all the others get further corrections 
from the integral term in the NLIE. Therefore for a general value of p the series ( [4.6| ) must 
be truncated to its first term for consistency. Using the asymptotic behaviour of Ki{z) we 
obtain 



M V Txl 



e~'cos79 + . . . . (4.^ 
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ln+l> 



T 



ln> 



X 



Figure 2: The one-instanton contribution to the vacuum energy 

When the twist angle takes the allowed values ( p. 8]) {}m = ^irm/k, this result can be 
compared to (|4.1|) , showing that the leading n = 1 term in (|4.6|) predicts the value of 
EiiL)/M. 

It is easy to interpret this result in terms of the usual instanton calculus. The one- 
instanton configuration is nothing else but the static one-soliton configuration of sine- 
Gordon theory on a cylinder replacing the space variable by the Euclidean time. It is 
independent of the spatial coordinate and so satisfies the periodic boundary condition on 
the cylinder. The Euclidean action is just the soliton mass M multiplied by the volume L. 
This gives us the factor e~' by the normal rules of instanton calculus. The factor ^j2/{'nl) is 
composed of two parts: a contribution of ^/^/{2^^) comes from the one (bosonic) zero mode 
of the instanton generated by translations and the rest from the determinant of the nonzero 
mode oscillations around the soliton, truncated to quadratic terms in the action. One then 
gets the result (|4.8|) from the usual dilute instanton gas calculation (cf. ( |A.21J )). (Since 
this determinant is different from the one needed to compute the quantum corrections to 
the SG soliton's mass we spell out the details of the calculation in appendix |^). Thus 
for Ei{L)/M the n = 1 term of the NLIE and the DIG give identical predictions. This is 
interesting, as the former one is expected to give reliable results for p ~ 1, while of the 
latter we expect this for p ~ 0; however, as we have seen the leading term is independent 
of p. Note that this prediction is also independent of the folding number k. 

To obtain a theoretical prediction for Ej{L)^ j > 1 in (|4.1| ), one must go beyond these 
approximations. In the NLIE — as mentioned earlier — this would necessitate the inclusion 
of the integral term, while in the instanton calculus it would require a handle on the exact 
multi-instanton solutions (as opposed to the approximate ones in DIG) together with their 
determinants in the cylindrical spacetime on Fig. |2|. To derive an explicit expression for 
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Figure 3: TCSA data versus predictions for E'l, E2 and -^3 at p = 2/7 



Ej{L) is beyond the scope of the present paper, and we merely note that we expect 



E,iL) 



C{3.P.k)-^Al + 0{e~^^^'^^)). (4.9) 



M > lm{j,p,k) 

Keeping simply the first few terms with higher n in (O) predicts 



c^o-,p,/^)=(-i)y^, 

and m = 1/2 for j = 2, 3, ... , and on Fig. |^ and Fig. ^ we used these values for mak- 
ing predictions about E2 and E^. Since the NLIE is an exact description, the constants 
C{j,p, k) and the exponents m can in principle be determined exactly (or numerically with 
very high accuracy). 

To test the NLIE and instanton predictions we determined numerically Ei, E2 and E3 
from the first different [k/2] TCSA eigenvalues using eqn. ( [4.1|) . On Fig. |^ we collected 



these quantities for models having the same p = 2/7 but differing in their folding number, 
which varied between 2 and 6. The data show a universal behaviour with no folding 
number dependence and in case of Ei they fit very well to the NLIE/instanton prediction, 
eqn. (|4.8|). Please note that the data and the most naive NLIE prediction ( |4.9| ) for E2 and 
Es differ only in the prefactor C{j,p, k) as the predictions run parallel to the data in the 
semilogarithmic plot. On Fig. ^ we compile the numerical values of Ei and E2 in 4- folded 
models having different values of p. Again the data show a universal behaviour with no 
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Figure 4: The p (in) dependence of Ei and E2 in various 4-folded models 

significant p dependence, and the prediction ([4 .81) for Ei describes the data very well. This 
figure also shows that the smaller p is the larger is the / range where TCSA gives reliable 
data. 

5 Multi-particle energy levels in finite volume 

The finite volume spectra of completely integrable models can be used to test their conjec- 
tured exact S'-matrices. In particular, the energy levels of multi-kink states satisfying the 
periodic boundary conditions can be determined in terms of their S'-matrices ||^. These 
energy levels then can be compared to the finite volume spectrum obtained by TCSA. 

There are three physical effects that contribute to the finite volume energy levels of a 
QFT. The 'tunnelling' effects are there in any theory — like SG(/5, k) — which has degenerate 
vacua in infinite volume. The corrections due to tunnelling are C(exp(— ML)), where L is 
the volume of the compact coordinate and M is a characteristic mass, in our case the mass of 
the quantum kink. In any massive theory there are two types of 'off-shell' effects due to the 
vacuum polarization and the interactions mediated by virtual particles, but both of them 
give C(exp(— ML)) corrections only. Finally in finite volume the particles continuously 
scatter on each other in a multi-particle state and as a result of these 'scattering' effects 
the stationary scattering states (i.e. the ones invariant under the mutual scattering of 
the constituents) are the true energy eigenstates. The resulting quantization conditions — 
called Bethe-Yang equations in @| — can be expressed in terms of the S'-matrix if we assume 
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that the particles are point-like. The corresponding corrections to the energy levels are 
0{L~'^) and so much larger than in the previous two cases, therefore they provide an ideal 
possibility to check the S'-matrices. 



5.1 Particle spectrum in the classical SG(/5, /c) in infinite volume 

As mentioned above the static soliton solution of the SG theory becomes a kink in the 
fc-folded model. More precisely the classical kink solutions of SG(/3, k), connecting neigh- 
bouring minima can be written with the aid of the SG soliton solution as 

ir„,„+i(x,t) = -arctan e^^^^^-^^o) + —- n = 0, . . . , k - 1; k = . (5.1) 

Thus instead of the single soliton we have k different 'one kink' solutions in SG(/5, k). (The 
antikink solutions are obtained by the (p t— > —(p reflection). Since a multi-kink solution 
corresponds to a sequence of vacua on a line one cannot arbitrarily compose single kinks 
to obtain an allowed solution, in contrast to the SG (anti)solitons. This restriction on 
the sequence of kinks translates into restrictions on the multi-particle Hilbert space in 
the quantized model. Note however that any multi-soliton/antisoliton solution of SG has 
a (multi-)kink interpretation in SG(/3, fc), at least in inflnite volume with no boundary 
conditions prescribed. 

SG(/5, k) also has k different breather solutions which oscillate around the k different 
minima of the potential: 

4 sin (^) 2n7T 

5^^) (x, t) = -arctan \^ +^ / + , veR, n = 0,...,k-l . (5.2) 

^ ^cosh(^) /5 

An important characteristic of the SG solutions is their topological charge. It measures 
how many times the (p{t,x) field winds around its range as x runs from —oo to oo. Thus 
in SG{P, k) we define it as 



Qk = ^ J d^ifdx . (5.3) 

This implies that the single kink solutions have a fractional 1/k topological charge and to 
make Qk integer we have to consider at least a /c-kink solution. The fractionally charged 
configurations give rise to nonlocal states at the quantum level, which explains the re- 
striction ( p. 91 ) imposed on the NLIE sources. The quasi-periodic boundary condition ( p7^ ) 
excludes all the single kink solutions, but the breathers ( ^.2| ) satisfy it, as do the multi-kink 
ones with integer topological charge (|0|) , at least approximately for L ^dass where 



A^ciass = is the classical kink mass. 



5.2 The particle spectrum of the quantum SG{/3,k) 

In infinite volume, as a result of its integrability, the quantized SG model contains particles 
corresponding to classical soliton or breather solutions. Since integrability is a consequence 
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of the local properties of the Lagrangian as a function of if we expect that the quantized 
SG(/?, k) also has kink and breather particles. Of course the symmetry requires that in- 
stead of the single quantum soliton (antisoliton) of SG in the fc-folded model there should be 
k quantum kinks (antikinks) degenerate in mass. In addition we expect breather particles 



-i-'n 



(m) 



n 



0, 1, ... /c — 1 corresponding to all the k different vacua. Furthermore, contem- 



plating e.g. a semiclassical analysis, we expect that the relation between the quantum kink 
mass M and the possible breather masses Mm is given by the familiar expression 



Mm = 2M sin 



TTpm 



m 



(5.4) 



independently of the 'vacuum' index n of bIi^\ This expectation is confirmed by Fig. |^ 
where, in the 2 folded model, for p = 2/7, we compare the breather masses (horizontal 
lines) with the TCSA data.0 These data are obtained with vanishing total momentum 
and vanishing topological charge, and the various dots represent the first eight energy 
eigenvalues above the ground state. The data corresponding to single particle states (i.e. to 
the first, second and third breathers), tend much faster to their infinite volume values than 
the two-particle lines, having 2Mi as their asymptotics. The reason is that while in the 
single particle masses there are only 0(e~') finite size corrections, in the energy of two- 
particle states there are 0{l~^) corrections coming from the mutual scattering among the 
particles. 

Denoting the kink of rapidity 9, interpolating between the vacuum a at x ^ — oo and 
/? at a; — oo by Kaj^iO), we let the amplitude 



describe the process 



07 \ 



+ K. 



Ka5 



where 9i > 62 and 612 = 61 — 62. Labelling the vacua by a = 0, 1 . . . , — 1 the model 
SG(/9, k) has kinks K^js with |a — /3| = 1 (or A; — 1). It is easy to describe the S-matrices of 
these kinks 0|: using the symmetry as well as time reversal and parity invariance one 
can show that every nonvanishing kink-kink amplitude is equal to one of the following three 
amplitudes: S'^^lT^^i^), ^"^^"+^(^12), -^aa ^""^^12), where a ± 1, a + 2 are understood 
mod k. These (a independent) amplitudes are also independent of the global properties 
of if and thus they should be equal to the soliton-soliton {Sss{6)), antisoliton-soliton re- 
flection {Sr{6))^ and antisoliton-soliton transmission {St{6)) amplitudes, respectively, of 
the ordinary SG model. To compute multi-kink energy levels we need the explicit form of 
Sss{6): 

du sm{6uj) sinh {\uj{p — 1)) 



Sss{6) 



m 



UJ 



cosh (f cj) sinh {^ujp) 

''On Fig. 5-8 E/M on the vertical axis stands for [E{1) — Eq{1)]/M, where Eq{1) is the vacuum energy, 
i.e. the ground state energy in the Q = sector. 
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Figure 5: Breather masses in the 2-folded model with p = 2/7 

In a completely analogous way one obtains that the S'ir\0i2) S'-matrix, which describes 
the possible breather-kink scatterings 

+ ^n,„,+l(^2) ^ ^n,„+l(^2) + SS(^l) , ^1 > ^2 



coincides with the breather-soliton S'i^V^i-?) S'-matrix of the SG model Wj 



oim, ... _ ^inh g + ^ cos ^ sin^ {^^p - | + zf ) 
^ ^ ^ sinh^-^cos^ sin2(ni^7rp-2_,|) ■ ^^"^^ 



In finite volume, when the boundary condition (|2.3|) is imposed the single kink particles 
disappear from the theory and they survive only as the building block constituents of the 
multi-kink state with integer topological charge. On the other hand the breather particles 
bIT'^ are there even in finite volume, as they are consistent with (|2.3|) . 



5.3 The Bethe-Yang equations 

To make a comparison with the TCSA results we need the multi-particle energy levels as 
functions of L. While it is possible to get them directly in the NLIE formalism, we shall 
use a simpler (and approximate) method which can be shown to give equivalent results for 
large volume (ML 3> 1) but which has a clear interpretation in terms of the spectrum of 
the model as it was described above. 
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In an integrable theory particle number is conserved, thus the concept of particle 
state with any fixed is well defined, at least for L ^ M^^ when most of the time 
the particles are far from each other. First we consider the 'pure' multi-kink states and 
envisage the iV-kink energy eigenstate in a large volume as a stationary scattering state, 
characterized by the set of (conserved) rapidities 9 = {9i, . . . ,9n}- Any A^-kink stationary 
wavefunction can be expanded using the independent allowed A^-kink m-states (the number 
of which we denote by (1^) 



\Kn,nM)Kn,ns{02) . . . if„^„, (^^7v) ) , 9, > 92 > . . . > 9 



N 



as basis, and let ip"-{9) (n = (ni . . . un)) denote its components with respect to this basis. 
Then, as a consequence of the periodic boundary conditions ip^{9) must satisfy the following 
Bethe-Yang equations (at least if all particles have the same mass) : 

giAfLsinhe, ^y.^^^m^n^^) _ _r^rn0^ ^ J = 1, . . . , A^ , (5.6) 
n 

where Tj{9) is the A^ particle transfer matrix 

TV 



=1 



For a given L eqns. (q^) have solutions only for some special 9j, and the total energy and 
momentum of the system in a state characterized by these solutions are given by 

N N 

E = ^Mcosh9j , P = J2Msinh9j. (5.7) 

i=i j=i 

Note however, that while this expression for P is exact even for finite L, the one for E is 
only approximate, as we neglect the 'tunnelling' and 'off-shell' corrections. 

It is straightforward to use this formalism to obtain the energy levels of pure multi- 
kink states in the Q = n {n > 1) sector of SG(/5, fc). It is natural to assume, that the 
lowest energy levels in this sector correspond to states with the smallest possible number 
of particles. The states with lowest number of particles compatible with the boundary 
condition and the topological charge being n contain N = k ■ n kinks and no breathers 
or antikinks. Since the sequence of kinks in these states is necessarily fixed, we have, 
independently of n, only k different basis vectors in this subspace, i.e. = k. The basis 
vectors can be chosen as 

\^o,i,...,k-i {9,,..., 9k.n)) = \Koi i9^) . . . K^-i,, (4) i^oi {Ok+i) ■ ■ ■ K^^ifi (Ok-n)) 
\i^k-i,o,...,k-2{9i, . . . ,9k.n)) = \Kk-ifl (92) Koi {9i) . . . 

Kk-2,k-l i9k) (dk+l) ■ ■ ■ Kk-2,k-l {9k.n)) 

|V'i,...,fc-i,o (^1, . . . , 9k.n)) = \Ki,2 {9i) . . . Ko,i {9k) i^i,2 {9k+i) ■ ■ ■ Ko,i {9k.n)) , (5.8) 
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where 6i > 62 > ■ ■ ■ > Ok-n- Using the explicit form of the transfer matrix and the 
identification of the various 5*2^ elements with the SG soliton-soliton S'-matrix Sss{0)-, the 
Bethe-Yang equations can finally be written in the following form: 



k-n 



AML sinh 6», 



1 k-n. 



where t/^ is a column vector made of the coefficients of the basis vectors in eqn. 
and the k x k matrix Tk 

/O ... 1\ 

1 ... 

n= 1 ... 



VO ... 1 oy 

describes a cyclic permutation generating Z^. The eigenvectors ipm belonging to the various 



eigenvalues (e*T 



i — m/k 



m 



0, . . . , /c — 1) of Tfc carry the different inequivalent irreducible rep- 



resentations of Zfc. Choosing the m-th eigenvalue of and introducing the dimensionless 
variable / = ML gives rise to the following Bethe-Yang equations 



1, . . .k ■ n , 



(5.9) 



where iVj G Z + 1/2 for ■ n even and Ni E Z for k ■ n odd, and as a result of Sss{0) 
we must have Ni 7^ A'^ for i j. The total momentum carried by a solution of (|5T 
determined solely in terms of the iVj: 



-1 
is 



P 

M 



271 

T 



m ■ n 



The simplest possibility — and the one we investigated by TCSA — is when the total (CM) 
momentum vanishes, P = 0. Please note that as a consequence of 5{—0) = —S{9) 
and e*^'^^''"'")/'^ = e~^'^™^'' the Bethe-Yang equations guarantee that the P = states 
{ipmyi^k-m), which Carry complex conjugate representations of D^, are degenerate in en- 
ergy. 

It is also possible to arrive at this result directly from the NLIE discussed in Section ||. 
In the infrared limit ML ^ 1 the integral term in the NLIE (|3.1|) becomes negligible and 
for a state with k ■ n holes the equation simplifies to 

k-n 



Z{9) =d + l sinh 9 + ^x{d - % 



Observing that x(6') = 8{9) and remembering that the allowed values of ^ (t3.8|) are exactly 
equivalent to selecting 



U r. 



m 

2.-, k 



0, 



m — 1 
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the Bethe quantization rules reduce to eqn. (^.91 ) while the energy/momentum formulas 
turn into (|5.7] ). Note that the quantization rule (p.7|) which selects the local operator 
algebra Ab from ( |2.19|) is exactly the one observed above for the iVj: for k-n even it assigns 



half-integer, while for k-n odd integer Bethe quantum numbers /j. We remark here that in 
general the large volume limit of the NLIE coincides with the Bethe- Yang equations (^) 
in their scalar form, i.e. evaluated on the eigenvectors of the transfer matrix T™(^). 

5.3.1 Comparison with the TCSA data 

To compare with the TCSA results we consider a few sectors of the models with k = 2 and 
k = 3 in more detail. The simplest of them is the Q = 1 sector of SG(/5, 2). In this case 
ipo describes the symmetric and ipi the antisymmetric wave functions of Kqi and Kiq, the 
latter one being allowed as the two kinks are different (bosonic) particles. P = implies 
in this case that 6i = —62 = and ( |5^ simplify to 



(.inh9 + *(20)=2.{;^»|^+|-_ (5.10) 

From this equation 6 = 6{l, N^) can be determined using e.g. an iterative procedure, or 
alternatively the volume dependence of the 2-kink energy levels can be given in parametric 
form cis 

This makes it clear that on the {l,E/M) plane the 2-particle lines cannot intersect each 
other for < / < 00. These findings make it possible to distinguish clearly between the 
solitons of sine-Gordon theory and the kinks in the 2-folded model, i.e. to argue that 
the Q = 2 sector of SG(/3, 1) (denoted as SG{P, 1)2) is different from the Q = 1 sector 
of SG(/3,2) (SG{P, 2)1). Indeed it is straightforward to derive the Bethe- Yang equations 
for the pure 2-soliton states in SG(/3, 1)2; since these solitons are identical particles these 
Bethe- Yang equations are given by the first line in eqn. (|5.10 ). Therefore the number 



of 2-soliton states in SG(/9, 1)2 is half the number of the 2-kink ones in SG(/3,2)i. On 
Fig. ^the continuous lines are given by the interpolated TCSA data obtained in SG(/5, 2)i 
with p = 2/7, while the dots represent the TCSA data obtained in SG(/9, 1)2; they clearly 
correspond to every second line onlyQ. 

On Fig. 1^ the continuous lines depending on one quantum number only are given by 
eqn. ( |5.10| ) and the dots now correspond to the TCSA data in SG(/5,2)i. On this figure 



we find data lines that cannot be interpreted as pure 2-kink states, as they apparently do 
intersect some of the other lines. Since the large I behaviour of these lines is compatible 
with E — > 2M + Ml for / 00, it is natural to try to interpret them as describing 
3-particle states containing one 'first' (m = 1) breather in addition to the two kinks. 
Since Mi 7^ M the transfer matrix formalism worked out in does not apply directly. 

^This is also true for the two lines which cannot be interpreted as two-particle ones due to their 
(multiple) intersections with the other levels. 
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Nevertheless following the original line of thought, namely by 'commuting around' any of 
the particles using the appropriate S'-matrices one can derive the Bethe-Yang equations 
for the stationary scattering states. 

Deleting the sequential or 'upper' (m) index of the breathers (since we consider only 
the first one) but keeping their 'vacuum' index the basis vectors in this 3-particle subspace 
can be chosen as 



V^o^i^ = 1^01(^1)^^10(^2)50(^3)) 

plus 4 other states corresponding to the breathers being 'in the middle' or in the 'front'. 
(The index h signals that the breathers are the last). Introducing the notation Ss^\6) = 
_Qia{e) ^ MiL = ml and = :^(^°^^ ± ip^^^)^ one can convert the Bethe-Yang equations 
derived by the 'commuting around' procedure into the following form: 

Ismh9i + 5{ei-e2) + a{6i-e^) = 2ttN^ (5.12) 
lsmhe2 + 6(62-61) + 0(62-63) = 2'kNI (5.13) 
mlsmh63 + a(6:i- 61) + 0(63-62) = 27rK^, (5.14) 

where 7 = ±, ^ ^1 ^12 ^ Z + 1/2, K'^ e Z. Once again, an equivalent system of 
equations can be derived from the NLIE in the ML ^ 1 (infrared) region. 

In the P = system there are solutions with the breather at rest and the two kinks 
moving in opposite directions: K'^ = = 63, N2 = —N^ 62 = —61. Then the (|5 .121 - 10^ ) 
system simplifies to the single equation 

/ sinh 61 + 6(26i) + a(6i) = 2ttN^, (5.15) 

which admits the parametric solution 

f27TN^ -6(26i) -a(6i) , , 

It is easy to understand how the = —N2 =1/2, K~ = state tends to the A^i = 1 
2-kink state (eqn. (|5.10| )) in the UV limit: in this limit a(6) — > — vr and this effectively 
converts the 1/2 on the right hand side of ( ^.15| ) into 1. The lines depending on three 
quantum numbers on Fig. ^correspond to the first two possibilities given by eqn. ( ^.15| ), and 
the agreement with the TCSA data is excellent. Note also that the line with A^^*" = = 1, 
= is also present in SG(/3, 1)2, consistently with the data on Fig. |. 
The Q = 1 sector of SG(/3, 3) is interesting, as the representations belonging to m = 1 
and m = 2 = — 1 are complex conjugate ones. The Bethe-Yang equations in this case take 
the form 

/ • ^ V" A n \ o / ^0'^ ^ ^ ' = ' h ^ = 1, 2, 3 and 

/smht^j+ > d [6i — 6j] = 271 < (■^ ^ where ^/i) , ^ . , . ■ 

j^, \ e Z ± i , m = ±1 , Nin' ^ N^Uoi t ^ J 
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Figure 8: Bethe-Yang curves and TCSA data for SG (/5 = |v^) 3) 



For the symmetric states (i.e. when m = 0) with P = there are special solutions with 
6*3 = — 6'i, 02 = (when Nq = —Nq = —N, Nq = 0); and for them the volume dependence 
of the energy levels can be given in parametric form similar to eqn. ( ^.llj ), but apart from 
these special cases we have to rely on numerical procedures to get -^(/,A^^). On Fig. ||, 
for p = 2/7, we show how well the predictions of these Bethe-Yang equations describe the 
TCSA data, regarding both the degeneracies and the volume dependence. 



6 Vacuum expectation values of local fields 

In this section we analyze the vacuum expectation values of exponential fields. Some time 
ago an explicit expression was given for this quantity by Lukyanov and Zamolodchikov [Q. 
Here we compare this prediction with the data extracted from TCSA, and by doing so we 
give further evidence that everything which in sine-Gordon theory follows only from the 
local properties of the scalar field and the Lagrangian remains true in the fc-folded model 
as well. 



6.1 Symmetries and vacuum expectation values 

In this subsection we will isolate the independent amplitudes which characterize the vacuum 
expectation values of exponential fields in SG(/3, k): 

{^n\ Vm \A) , Vm = V(mfl) = exp U(3 — l^] . 
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Our considerations will be valid for any value of the volume L. 

As a consequence of the action of T on the vacua (|2.10|) and on the local fields Vn-, 
( p.l4| ) we have 



(^„| Vm \^r) = A^{r)Sn,n.+r , (n| V^ \r) = B^{r - n)e^^^ , (6.1) 

where Am{r) are certain (unknown) amplitudes and 

fe-i 

Bm{r) = (0| Vm \r) = -^Yl e~'^"''A,n(n) . (6.2) 

n=0 

Reality of the field (p yields 

Vl = V.m 

and so 

BmivY = e¥-S_„(-r) . 
The transformation properties of the field ip under S imply 

Emir) = 5_^(-r) . 

This allows us to determine the phase of the amplitude Bmij) up to a sign: 

Bm{T) = ei'^™F^(r) , F™(r) G M . (6.3) 

It turns out that the real amplitudes Fmij-) are not all independent. Indeed, as it was 
already remarked in subsection |2]^, in the usual basis (|2.21j ) of the (ultraviolet) free bo- 
son Hilbert space all the Vm have real matrix elements. Therefore the Hamiltonian (^) 
as a matrix is real and symmetric, and as a consequence all its eigenvectors have real 
components. This implies 

mVm\A)eR 

and result 

Fm{r) = {-irFm{k-r), r = 1, . . . , A; - 1 . (6.4) 

Note that there is no constraint on -Fm(O) = (0| |0). The vacuum expectation values 
can therefore be characterized by the independent amplitudes 



Fm{r) , r = 0, . . . 



(6.5) 



all of which are real. 
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6.2 The Lukyanov-Zamolodchikov formula 

In § Lukyanov and Zamolodchikov proposed an exact formula for the vacuum expectation 
values of exponential fields in the ordinary sine-Gordon model in infinite volume. Here we 
briefly recall their result. We normalize the exponential fields so that the short distance 
asymptotics of the non-vanishing two-point functions reads 



(0|e 



iaip{x) -iaip{y) 



|0) 



\x — yP^" 



where = |^ is the conformal weight of e*""^*^'^-', and |0) denotes the state \n) for n 
(see eqn. (|275|) ). Let us define 

g{a) = (0|e'"'^(") |0) . 

The authors of § conjecture^ 



g{a) 



exp 




si 

47r 



sinh^ ifj) 



-2t 



which is valid for 



_2sinh(^t) sinh(t) cosh 



< Stt and |3fte a(3\ < Att , 



An 



(6.6) 



(6.7) 



and where M is the soliton (kink) mass. Recalling now our basic idea, namely that every- 
thing which depends only on the local properties of and the SG Lagrangian also holds 
for the fc-folded model SG(/?, k), we can identify Q{a) and -Fm(O): 



m(3 



F„(0) when I = ML = oo. 



The TCSA data allow us to extract at finite I the expectation value of e*'"^*^"'''"^||^|=i (where 
w = is the map from the cylinder to the plane) in the ground states found by the 
numerical diagonalization, i.e. we can measure Am{n) as functions of /, 74m(^^)[/]. Thus 
introducing the dimensionless function g{a) = M^^^°-Q{a) we obtain for finite /: 



n=0 



^^g{a)N{l), a = ^ 



(2vr) 



(6.8) 



where A^(/) is a finite size "correction" factor of which we know N{1) ^ 1 + C(e ') for 
/ oo. 



'Our normalization for the field (p and the coupling constant f3 differs from that of S] by a factor of 
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The amplitudes Fm{r) are in general not known for r 7^ 0. However, one can derive 
their leading behaviour for large volume using the following consideration. First note that 
the matrix elements 

{0\Vm\r) 

vanish in infinite volume since there is a super-selection rule making the vacua \n) lie in 
physically disconnected Hilbert spaces with no local operators connecting them. It is clear 
that in finite volume the non-vanishing contribution to these matrix elements comes from 
the vacuum tunnelling described by instanton effects, whose magnitude we calculated for 
large I in section ^. Therefore we expect 

(0| Vm \r) ~ e"' when / > 1 and r 7^ . 



6.3 Comparison to TCSA 

We tested the Lukyanov-Zamolodchikov formula in a 3-folded model with p = 2/7. Using 
the eigenvectors belonging to the three ground states found by the TCSA algorithm we 
determined numerically the various amplitudes for all the values of m (m = 

1, . . . , 5) satisfying the conditions (|6.7|) . Then, recalling eqn. (|6.8| ), we plotted Fm(0)[/]/~^^™ 
where 

2A = — 

as a function of 1. The result is shown on Fig. |, where the horizontal lines correspond to 



m(3 



The agreement between the measured values and the predicted ones is very good, though 
for higher values of m the data are somewhat below the horizontal lines indicating that 
the TCSA data should be extrapolated as in the discussion that follows. The 3-folded 
model provides a good laboratory as it makes possible to extract from the data 

as well. The result Fi(l)[/] can be seen on a semilogarithmic plot on Fig. |10|. Clearly this 
behaviour is consistent with the expected exponential fall-off. 

The exact expression for Q{a) was obtained in by a clever interpolation between 
various limits, where the corresponding expressions were known from other sources. One 
important such limit was the semiclassical one. Since our numerical study is for p < 1, 
we investigated whether using our data one can make a distinction between the exact 
and semiclassical expressions of i.e. whether one can justify the exact or merely the 
semiclassical formula. 

We carried this out by zooming in on the vicinity of the uppermost (m = 1) line 



on Fig. y, and the results are compiled on Fig. |TT|. Here the upper/lower horizontal 
line corresponds to a g{a,l3) obtained by the exact /semiclassical expressions in [|], and 
the TCSA data displayed were taken at four different -Ecut values. At each fixed / the 
data apparently converge monotonously with increasing -Ecut- The validity of the exact 
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Figure 9: TCSA data and the Lukyanov-Zamolodchikov prediction 
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Figure 11: Comparing TCSA with the semiclassical and full quantum formulae 



expression over the semiclassical one is supported by the fact that these monotonously 
convergent data exceed the semiclassical expression in a certain range, while they always 
stay below the exact one. To strengthen this conclusion, at each fixed /, we extrapolated 
the data by fitting their i^cut dependence with an expression 

,n gS'"'+c(i) 

and the extrapolated points on Fig. |ll] correspond to the coefficients a(/).[] These extrap- 
olated points show a monotonously increasing behaviour in / and they significantly exceed 
the semiclassical expression while always stay below the exact one. In the final step, re- 
calling eqn. ( |6.8|) , we fitted the extrapolated points hy A — Be~^ /I, where e~' is motivated 
by the instanton contribution, and the by the fiuctuation determinant without a zero 
mode. This 'instanton fit' describes the extrapolated data with a very small variance, and 
though the A obtained this way is somewhat smaller than the exact expression, it is much 
closer to the exact than to the semiclassical one. To sum up, we can say that our data 
indeed favour the exact expression of over the semiclassical one. 

We remark that a similar calculation of vacuum expectation values from TCSA was 



performed in ||T8| for the case of $(1,3) perturbations of Virasoro minimal models. These 



^This extrapolating formula is consistent with the monotonous increase and the power of i?cut in it is 
motivated by the observation, that envisaging the determination of G{a) in pCFT the effective expansion 
parameter would be rather than A. 
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models are restrictions of sine-Gordon theory at rational values of p and as a consequence 
one can derive a formula for the vacuum expectation values of local fields starting from 
( |6.6| ) H). Similarly to our case, they find that the respective formula agrees with the TCSA 
results. In contrast to their approach, however, we checked the prediction ( |6.6|) directly in 
the case of sine-Gordon theory. 

There is another important implication of this result. The Lukyanov-Zamolodchikov 
conjecture for the vacuum expectation values of local fields is connected to (and can be 
derived from) the conjecture for the so-called Liouville refiection factor ||19[ |20|. The above 
verification of the formula ( |0|) therefore lends an indirect support to the conjectured 
expression for the refiection factor and considerations based on it. 



7 Conclusions 

In this paper we investigated the fc-folded sine-Gordon model SG(/3, k) in finite volume. 
The aim of this study is to give support to the idea that the fc-folded boundary conditions 
which make SG(/?, k) different from the ordinary sine-Gordon theory do indeed preserve 
the integrability of the model, while changing the spectrum in a well defined manner. 

We analyzed three major problems in some detail and showed in all of them, that the 
consequences one can draw from this expectation are indeed correct. In particular we found 
that the leading part of the split in the ground state energy levels, for which the NLIE 
and the instanton calculus gave identical results, does indeed coincide with the numerical 
data obtained by using TCSA. Furthermore, we provided evidence that the k dependent 
degeneracies and the volume dependence of the multi-particle energy levels can indeed be 
described by the formalism of |Q, thus indirectly we verified the conjectured S'-matrices. 
Last but not least we showed that the vacuum expectation value of the exponential field 
can be measured in the k folded model and we gave evidence supporting the validity of the 
formula proposed in 
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A Instanton calculus in finite volume 

In this appendix we perform a standard instanton calculation based on the dilute instanton 
gas approximation, following the outlines of ||2T]|. Here we shall work in a /c = cxd theory 
which means that we drop the identification of the field in (|2.2| ). We will comment on this 
issue later. 



A.l The dilute instanton gas approximation 

The Euclidean action of sine-Gordon theory is 

foo />L/2 „ 1 . //2 

I 

-L/2 

The equations of motion following from the action ( |A.1J ) admit the one-instanton solution 



Se[^] = dr dx Q (d^^f + \ {d^^f + ^ (1 - cos/3^)) . (A.l) 



4 

V'inst = arctanexp (/io (r - tq)) (A. 2) 

with action ^ 

Se binst] = = Mci^^L , (A. 3) 

where Mdass is the classical soliton mass. Using the usual rules of instanton calculus, the 
saddle-point evaluation of the Euclidean path integral yields the following result for the 
level splitting as a function of the twist angle t9 labelling the energy eigenstates 



^E{d) = -2cos(i9) 



det'M 


-1/2 


^ Se [v^inst] A 


detMo 


1 


V 27r J 



1/2 



where M and Mq are the operators 

M=-dl-dl + V" ((/.inst) , Mo = -df -dl + V" (0) (A.5) 

describing the fiuctuations around the instanton to quadratic order, det' M denotes the 
determinant of M without its zero mode and the factor 



Se[v, 



-"inst 



2n 



(A.6) 



comes from the one translational zero mode tq of the instanton solution ( |A.2| ). The cos^ 
dependence arises from Fourier transforming the dilute instanton gas summation according 
to the relation between the vacua |n) and \{}): 

oo 
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Here 1!} can take any value: the physically inequivalent choices are — vr < i9 < it. For 
SG(/3, fc), the result of the dilute instanton gas calculation is similar to (|A.4| ) with the only 
difference that "i? can only take the discrete values 



k 



mod 2tt . 



A. 2 The heat kernel representation for the determinant 

Let us define the heat kernel for a positive hermitian operator A in the following way: 

Kt{A)= Tie-'^, (A.8) 
from which the determinant can be reconstructed: 



log det A = Tr log A 







°° dt 



(A.9) 



As the determinant is divergent (the divergence comes from the lower end of the integration 
over t), we shall compute the difference 

Kt{M) = Kt{M)-Kt{Mo) 

using (^-function regularisation. We define the ^-function 

1 



Viz) 



dt t'-^ {Kt{M) - 1 



Then the determinant can be expressed as 



det'M 



detMn 



"1/2 



exp(lAc(.,M) 



(A.IO) 



(A.ll) 



z=0 



We can separate the x and r dependence and rewrite the heat kernel in the following form: 
kt{M) = Kt {-dl) {Kt (Q) - Kt (Qo)) , 

V cosh {ij,o{t-To)) 

Qo = -d'r + d- 
The heat kernels for the three operators that appear above are 

L — -2r2 



Kt {-dl) 
Kt{Q) 



E 

1 + 



^-n)t 



/Ant 



E 



n L 

e 4t 



dpp(p)e-*(^°+P' 
dppo{p)e-<^o+p') . 



(A.12) 



where p, po are the spectral densities for the operators Q, Qq and the additive 1 in Kt{Q) 
comes from the zero mode. 
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A. 3 Evaluation of the spectral densities 

The spectral density of the operator 



Q = -9r+f^l(^-^;^--) (A.13) 

V cosh (/Uot)/ 
can be evaluated by solving the spectral problem: 

- ^" - ililiM^ + fil^ = , (A. 14) 

cosh [fioT) 

where for our case s = 1. This problem can be solved exactly by mapping the above 
equation to a hypergeometric one p2| . 

A. 3.1 The discrete spectrum 



The discrete spectrum of (|A.14|) corresponds to < A < /Iq. The condition for square 
integrability reads 

1-^-5 = -n, nEN. (A. 15) 

For our case (s = 1) the only solution of ( [A.15| ) is n = which corresponds to A = 
i.e. exactly the unique zero mode of the operator Q mentioned before. (The other possibility 
n = 1 means A = /ig which is where the continuous spectrum starts.) 

A. 3. 2 The continuous spectrum 

The continuous spectrum covers the range fi^ < X < oo. In this domain we define a solution 
with the following asymptotic property: 

^(r^oo)=e^P^ X = fil+p\ 

Then one can compute 



^ ^ r(-.)r(i + .) r -^-s]t 1 + 

For s G N the first term vanishes, which means that the potential is refiectionless. For 
s = 1 we obtain 

v[;(r ^ -oo) = e^(^-+^(^')) , e^%) = !^±i^. (A.17) 

tp - /io 

To get a well-defined spectral density, we must put the system in a large box of size T. 
Then periodic boundary condition on \i/ implies 



pT - S{p) = 27iN , G Z. 
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For the size T — > oo the density of states with a given p is 



which gives 



pip) = TT ( ^ ~ 2^!^° 2 



A similar, but much simpler reasoning for Qq yields 

T 

Poip) = TT ' 

SO the divergent parts linear in T drop out from the difference. 



(A.18) 



(A.19) 



(A.20) 



A. 4 The instanton contribution 



Substituting the results ( |A.19| , [Ar20| ) into ( |A.12| ), the final result for the C-regularized heat 
kernel (|A.10D is 

Ciz, M) = Ci(z, M) + C2(^, M) + C3(^, M) , 

where 



Ci(z,M) 



—- / dte-'L{Antr'/^ {K,{Q) - K, (Qo) - 1^ 
J- {z> Jo 

d»-(A-.(-a^)-i), 



1 / \ 

^ dt t^-i (i^i - L (Ant)-'/') iUQ) - K, (go) - 1) 



Ci and correspond to renormalizing the classical soliton mass Mdass to the quantum one. 
The interesting contribution to the determinant comes from the term (2 which gives 



dC2{z,M) 



dz 



2=0 



Ant 



5 -lo£ 



An 



Numerical evaluation of the integral gives 



B 

2e~ 



TT 



to very high precision. Collecting all terms we get 

AE{^) = -M\ —e-^ cos^. 

V Til 



(A.21) 



REFERENCES 



35 



References 

[1] A. B. Zamolodchikov, Adv. Stud. Pure Math. 19 (1989) 641. 



[2] S. Ghoshal and A. Zamolodchikov, Int. J. Mod. Phys. A9 (1994) 3841, h^- 



th/9306002 



[3] J. A. Swieca Phys. Rev. D13 (1976) 312. 

[4] T. R. Klassen and E. Melzer, Nucl. Phys. B382 (1992) 441. 

[5] V. P. Yurov and Al. B. Zamolodchikov, Int. J. Mod. Phys. A6 (1991) 4557-4578. 

[6] A. Kliimper and P. A. Pearce, J. Stat. Phys. 64 (1991) 13; 

A. Kliimper, M. Batchelor and P. A. Pearce, J. Phys. A24 (1991) 3111. 

[7] C. Destri and H. J. De Vega, Phys. Rev. Lett. 69 (1992) 2313-2317. 

C. Destri and H. J. De Vega, Nucl. Phys. B438 (1995) 413-454, |hep-th/9407TT7 . 



[8] S. Lukyanov and A. B. Zamolodchikov, Nucl. Phys. B493 (1997) 571-587, |e£ 
th/9611238|. 



[9] T. R. Klassen and E. Melzer, Int. J. Mod. Phys. AS (1993) 4131-4174, |hep-th/9206TT4 . 
[10] G. Feverati, F. Ravanini and G. Takacs, Nucl. Phys. B540 (1999) 543-586, [hep 



th/9805117 . 

G. Feverati, F. Ravanini and G. Takacs, Phys. Lett. B444 (1998) 442-450, [lieF 
th/980716a 



[11] Al. B. Zamolodchikov, Int. J. Mod. Phys. AlO (1995) 1125-1150. 

[12] Al. B. Zamolodchikov, Nucl. Phys. B432 (1994) 427-456. 
Al. B. Zamolodchikov, Phys. Lett. B335 (1994) 436-443. 



[13] G. Feverati, F. Ravanini and G. Takacs, Nucl. Phys. B570 (2000) 615-643, [hi^ 



th/9909031 



[14] P. Fendley and H. Saleur: Massless Integrable Quantum Field Theories and Massless 
Scattering in (l + l)-Dimensions, [hep-th/9310058| . Lectures given at Summer School 
in High Energy Physics and Cosmology, Trieste, Italy, 1993. Published in Trieste HEP 
Cosmol. 1993, pp. 301-332. Also in Strings 1993, pp. 87-118. 

[15] D. Fioravanti, A. Mariottini, E. Quattrini and F. Ravanini, Phys. Lett. B390 (1997) 
243-251, [hep-th/ 9608091 . 



[16] C. Destri and H. J. De Vega, Nucl. Phys. B504 (1997) 621-664, [hep-th/9701107 . 
[17] A. B. Zamolodchikov and Al. B. Zamolodchikov, Ann. Phys 120 (1979) 253. 



REFERENCES 



36 



[18] R. Guida and N. Magnoli, Phys. Lett. B411 (1997) 127-133, |hep-th/97060T7| . 

[19] A. B. Zamolodchikov and Al. B. Zamolodchikov, Nud. Phys. B477 (1996) 577-605, 
[hep-th/950613q . 

[20] V. Fateev, S. Lukyanov, A. B. Zamolodchikov and Al. B. Zamolodchikov, Phys. Lett. 
B406 (1997) 83-88, |hep-th/9702f90[ . 

[21] G. Miinster, Nucl. Phys. B324 (1989) 630. 

[22] L. D. Landau and E. M. Lifshitz, Theoretical Physics lU: Quantum Mechanics, 
Butterworth-Heinemann, 1997. 



